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We present a new, simplified analysis of the low-energy electron-positron interaction, and use 
the resulting effective theory to calculate the binding effects that contribute to the decay rate, F, of 
orthopositronium through 0(c? lnaF). We express the total decay rate in terms of the annihilation 
rate for a free electron and positron at threshold, which has just recently been computed to sufficient 
precision. Our 0(a 2 F) result corrects errors in a previous analysis. 
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There has been longstanding uncertainty regarding a possible disagreement between theoretical and experimental 
determinations of the orthopositronium ground-state decay rate [jjj . A discrepancy here could have important implica- 
tions. For example, it might indicate a failure of perturbative expansions in a, the fine structure constant of quantum 
Cn| ! electrodynamics (QED); or it could signal the presence of new physical phenomena beyond QED. The decay rate is 
currently of particular interest because important parts of the theoretical calculation have recently been completed: 
in Ref. Q the annihilation rate for a free electron and positron at threshold is calculated through 0(ma 5 ), where m 
is the electron mass. This is two powers of a beyond leading order. Here we convert this decay rate for a free electron 
and positron into the decay rate for a positronium atom, the electron-positron bound state. We also calculate new 
, contributions at C(a 3 lnar). 

Traditionally, precision bound-state calculations have been formulated within the context of quantum field theory. 
, A Bethe-Salpeter analysis is an example; more recently, the nonrelativistic QED (NRQED) effective field theory 
approach has often been used It is, however, much simpler to recast the problem within the familiar framework of 
| nonrelativistic quantum mechanics || . Long-range QED interactions correspond to standard long-range potentials in 
the nonrelativistic Hamiltonian, while short-distance effects are described by a small number of local operators whose 
coefficients must be determined by comparing with the complete, relativistic formulation of QED. Here we relate the 
short-distance coefficients in the positronium Hamiltonian to the threshold annihilation rate computed in jlj , thereby 
obtaining the full Hamiltonian needed for computing the decay rate. The decay rate enters our formalism through 
the nonhermiticity of our Hamiltonian and the complex energies that result: T = — 2Im_E. 

The methods used in this paper are new to QED applications, and have a number of desirable features. A finite 
I ultraviolet cutoff is built into our Hamiltonian. This cutoff excludes high-momentum states that are poorly described 
by the nonrelativistic dynamics. Unlike traditional approaches we do not take the cutoff to infinity at the end of the 
O i' calculation; rather it is held fixed at a value of order the electron's mass. Consequently, no divergences occur, and the 
D . resulting energy eigenvalue problem can be solved nonperturbatively using simple numerical methods — for example, 
' by evaluating the matrix elements of the Hamiltonian using a finite basis set, and then solving a matrix eigenvalue 
, problem. This eliminates the need for bound-state perturbation theory. Another feature is that only physical inputs — 
on-shell scattering amplitudes, for example — are required from full QED, and therefore gauge and QED-regulator 
independence are explicit. Finally, our approach can be adapted in a natural way to multi-electron systems such as 
helium, where the lack of an exact zeroth-order solution and the complexity of bound state perturbation theory makes 
a nonperturbative solution particularly convenient. 
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I. THE EFFECTIVE HAMILTONIAN 



We now proceed to the effective Hamiltonian. We work in the center-of-momentum frame of the electron and 
positron, and consider only states of orbital angular momentum (5-states), and spin 1 (triplet- or ortho- states). 
In order to make use of the threshold results in JjJ, we give the photon a small mass, A, which is taken to after the 
determination of the local operator coefficients; the three final-state photons are not given a mass. 

We begin with Hamiltonian 

H^t---^ + V + iW (1) 

m 4m' 3 

where V and W are hermitian. Potential W accounts for the effects of three-photon decay, while V accounts for all 
other interactions. To determine corrections through 0{a 2 T), or 0(a 3 T), we need retain only non-annihilation terms 
that contribute to the positronium S -state binding energy through 0(a 4 m), or 0(a 5 m), respectively. We write V 
as the sum of three terms, 
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V(E) = V + V iel + V iad (E), 
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where the Coulomb potential, Vo, and the leading relativistic corrections, V Te \, are given by: 
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We have introduced Gaussian factors to suppress the potentials at high momentum. Typically we take the ultraviolet 
cutoff A ss to, although our final decay rate is (and must be) A- independent. With potentials Vo and V rc \ in place, our 
Hamiltonian correctly reproduces the QED amplitude for one-photon exchange and one-photon virtual annihilation 
to lowest and first order in (v/c) 2 (of the electron). V va< i(E) accounts for one-loop radiative corrections: 
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V[. ac j (E) gives in positronium the analogue of the non-recoil Lamb shift and the recoil Salpeter correction in hydrogen. 
The final term is a local operator accounting for effects at high momentum; parameter D is a counterterm which will 
be determined shortly. 

The annihilation of the electron and positron occurs over distances of order Ax ~ 1/to, which are much smaller 
than wavelengths typical of the electron and positron in the atom, A « 1/am. Thus the annihilation potential W 
consists entirely of short-distance interactions, which, to the order of interest, can be parameterized for 5*-states as 
follows: 



(i|W|fc) = A<® 
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We will adjust parameters A^°\ A^\ A^ and so that our Hamiltonian reproduces QED results for electron- 
positron annihilation into three photons. Since determination of A^ 1 and of B^> requires the as-yet unknown 0(ma 6 ) 
threshold annihilation rate, and the leading term in the momentum expansion for the 0(ma 4 ) rate, respectively, we 
simply set A^ = B^ 1 = 0. Doing so introduces an error in the decay rate of 0(a 3 T). 

The p 4 operator in Eq. ([!]) is ill-defined at high momenta. To regulate this operator, we replace it by an energy 
dependent potential ||: 
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Keeping only the relevant parts of this expression, our final effective Hamiltonian is therefore: 

H(E) = H (E)+iW(E) 



where 
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Potentials V (Eqs.(|),(|),(|),(|)), V (Eq.@) and W (Eq.(|)) are given above, and 
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(Z|Wo|fe> = (l\Wi\k)/A^ = A^eSST. 

Parameter D in Eq.(^|) is tuned to correctly reproduce the one-loop contribution to 
and is found to be 
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The hermitian Hamiltonian Hq(E) is accurate through 0(a 5 m). 

The parameters in W are determined by considering the imaginary part of the iS-wave scattering amplitude for 
ee — ■+ 37 — ► ee, with electron momentum k in the center-of-momcntum frame; the optical theorem relates this 
amplitude to the free-particle annihilation rate. For small k, the imaginary part of this amplitude can be parameterized 
as: 
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Parameters 7q, /?, ao, ai, 60 are 



We define T(fc) using nonrelativistic normalization for the external particles j 
determined using QED perturbation theory ||. 

We now calculate T in our Hamiltonian theory, and adjust the unknown parameters to reproduce the QED result 
order by order in a and k 2 . Matching at lowest order in a implies that 
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and A^> are determined by matching the 0{a) and 0(a 2 ) contributions in T(0), respectively; we obtain 
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Here 7 
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0.577216 is the Euler constant and k = j dx \nx 1 exp (— x 2 )erf (x) 2 — 0.051428. Having 



determined all the necessary parameters, we can now safely set A = 0. 



II. THE DECAY RATE 



Now that our Hamiltonian is completely specified we finally solve for the decay rate, given by the imaginary part of 
the ground state energy eigenvalue. Note that due to the presence of the cutoff, no divergences occur when calculating 
matrix elements, and no intricate limiting procedures are necessary to solve the eigenvalue problem — renormalization 
is automatic. To avoid dealing with nonhcrmitian matrices we choose to work only to first order in the annihilation 
potential W\ higher-order terms are suppressed by several powers of a beyond the precision of interest. 

We first solve the eigenvalue problem for Hq, 

H (E )\iPo)=E \il>o), (17) 

to obtain the ground-state energy and wavefunction. The energy dependence of Hq is easily handled by iterating the 
eigenvalue equation, starting with an approximate energy in Ho; the answer converges to adequate precision after 
only a few iterations. The eigenfunctions for our energy dependent Hamiltonian must be normalized so that [g 

= 1. (18) 
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Then the decay rate, to first order in W, is 
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rate. The cutoff operators v\(r) and S\(r) = — V 2 v\(r)/A-K are defined by Fourier transform 
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The matrix elements M\ , M 2 , in Eq. ( p0| ) can be evaluated for any S-state of H and the corresponding decay rate 
computed using this equation. We used bases consisting of 30 to 60 gaussians, with varying widths, to diagonalize H - 
The numerical eigenvalues accurately reproduce the 3 Si spectrum through 0(a 5 m) ||. Our ground-state results, for 
several values of A, are shown in Table |; there we introduce the dimensionless parameter Ar, defined for any S'-state 
by: 



T(nS) = ^ 
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This definition anticipates the leading a 2 In a Q and a 3 In 2 a |T(J contributions, which are correctly reproduced by 
our numerical analysis. The final results for Ar are almost independent of A, while the changes in the matrix elements 
from one A to the next are two orders of magnitude larger than Ar itself. Renormalization theory guarantees that 
A dependence due to the matrix elements cancels, in the final answer, against A dependence due to coefficients A^\ 
A^ and in the annihilation potential W. The residual A dependence in Ar is due to our approximations in 
potentials V and W, the dominant effect being at 0(a 3 T) where we have left out the contributions from ^4^ 3 ^ and 
£?«. 

Our calculation is nonperturbative in potential V and so automatically includes order a 3 T (and higher-order) 
corrections to the decay rate. To facilitate comparison with other calculations, we suppressed these higher-order 
effects by extrapolating our analysis to a = 0. Our results for A = m and a range of a's are shown in Table O. Upon 
extrapolation to a = 0, we obtain: 

( 0.913, n= LS 

x r(n)\ a ^ Q = < 2.588, n = 2S (23) 
[ 2.936, n= 3S 

where all results are accurate to within 1 in the last digit. The IS result agrees well with the analytic result in Q). 

Since the analysis would be complete through 0{a 3 T) with the inclusion of the operators parameterized by A^ 
and and since neither of these operators generates factors of lna, the complete contribution in the decay rate 

of 0(a 3 lnaT) is already present, along with that of 0(a 3 In 2 aT). By examining the a dependence of Ar, we find: 

Ar(lS) = 0.913- 0.665 a lna + 0(a). (24) 

The coefficient of a In a is independent of the cutoff A. The non-logarithmic term of 0(a 3 T) is cutoff-dependent, since 
we have neglected contributions from the cutoff-dependent and Using the values for ao and Tq in Ref. [0, 

the C(a 3 lnaT) contribution amounts to 2.4 x 10 _5 /is _1 . Together with a small contribution from five-photon decay 
of 0.73 x 10 _5 /is _1 |ll|], this brings the current theoretical prediction for the decay rate to T = 7.039967(10)/is~ 1 . 

After completing our calculation, we learned of an independent analysis of 0(a 3 mal 1 ) contributions to positronium 
decay by Kniehl and Penin Jrlfl . Our result, —0.665 for the coefficient independent of ao, disagrees with their 
result, —(4/5 + 81n2/3)/7r = —0.8430. To verify our analysis, we compared it with published results on C(a 3 lna) 
contributions to muonium hyperfine splitting (HFS). The HFS results involve the same operators, with different 
coefficients, as those contributing at the same relative order in the decay rate. Terms in the orthopositronium decay 
rate not involving ao arise from second order perturbations of V Ia d with the leading decay operator: 

ST jV rad GS 3 (r)) dV Iad 

r = 2 (S 3 (r)) +( ^ } ' (25) 

where G is the Coulomb Green's function with the ground state pole removed, and expectation values are taken 
between unperturbed Coulomb eigenfunctions. V Ta ,d can be expressed as 
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The reduced mass m r equals to/2 in positronium. The logarithmic parts of the matrix elements for the second-order 
perturbations can be inferred, for the most part, directly from the HFS papers [|l3| |l4| ]: 

-41n 2 a- 8(-ln2 + 3/4)hm , i = l 

ln 2 a+(21n2-l)lna ,i = 2 (28) 
2 In a , i = 3 

The only exception is the coefficient (2 In 2 — 1) of In a for C 2 ■ A partial analysis in Ref. Ill] gives (2 In 2 + 2) in place 
of (2 In 2 — 1). We have calculated the full contribution analytically, and find the result shown in Eq. ( ps| ) jlH]. We 
have also verified the results in Eq.(p8[) by direct numerical evaluation, using our Gaussian basis set. Non-logarithmic 
terms of O(a , a 1 ) are cutoff dependent, but the lo gar ithmic terms were cutoff independent, as expected. Combining 
our analytic results with those in Ref. fjj], our Eq. ( p4f ) becomes: 

X r (lS) = H _ 2 ln 2 + ^ 8 ln 2 _ ^ 2 In a (29) 

= 0.9129- 0.6647a ln a, 

in complete agreement with our numerical analysis. 

We thank the authors of [jjj for sharing their results with us before they were published. We also thank Patrick 
Labelle for several discussions. This work was supported by a grant from the National Science Foundation. 



A/to Mi M 2 X r 

0.5 18249.0812 37.6279 0.8518 
1.0 18405.3892 75.8182 0.8867 
1.5 18410.4321 113.7030 0.8323 
2.0 18375.7156 151.2649 0.5852 



TABLE I. Matrix elements and corrections to the ground state decay rate evaluated at a 1 = 137.03599976 as a function of 
cutoff A. To determine Xr, we set ao — a\ = bo = 0. 



a 


Ar(lS) 


X r (2S) 


A r (3S) 


0.08 


0.5362 


2.1012 


2.4437 


0.04 


0.7299 


2.3607 


2.7094 


0.02 


0.8282 


2.4870 


2.8373 


0.01 


0.8749 


2.5449 


2.8949 


0.005 


0.8962 


2.5703 


2.9198 


0.0025 0.9057 


2.5812 


2.9305 


-> 


0.913 


2.588 


2.936 



TABLE II. Corrections to decay rate at A = m as a function of a. Here A m = a + 2.444821528 
and A (2) = b - 2ai + 2.444822a + 6.179923. 
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